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Trigonal distortion of topologically confined channels in bilayer Graphene. 
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In this work we show that the trigonal warping of the electronic bands in bilayer graphene dramatically 
modifies the behavior of the topologically confined one-dimensional modes due to an inhomogeneous bias 
that changes sign across a channel. Up to four zero modes are present, depending on the orientation of the 
channel, these zero modes are related with a fractionalization of the topological charge due to the trigonal 
warping. 



Intervalley transport and coherence in graphene-based 
systems is an appealing subject, due to both the in- 
triguing physical phenomena that it entails and the 
great potential that it displays from an applied per- 
spective. Among others, valley-dependent resistance ef- 
fects and valley-filter devices, have been proposed in sev- 
eral contextsii^. Importantly, basic theoretical estimates 
prove that such devices are well within current or soon 
to come technological capabilities. 

Bilayer graphene (BLG), two coupled graphene lay- 
ers, has an electronic structure remarkably different from 
graphene. Unlike single layer graphene (SLG) the elec- 
tronic dispersion near the K points is parabolic and a 
gap can be induced by applying an external electric fieldS, 
perpendicular to the layers. 

The coherence between carriers on each of the layers of 
BLG opens up an additional effect. It has been shown^ 
that biased BLG is characterized by non-trivial topologi- 
cal structure associated with the behavior of the pseudo- 
spin degree of freedom. This topological property leads 
to several fascinating predictions j^"— In particular it can 
be shown that confined gapless bands can be induced in 
BLG by applying an inhomogeneous bias that changes 
sign across a channel. In this way a one dimensional 
channel populated with chiral zero modes is created with 
potential applications. It can be done for instance in a 
graphene p-n junction^, provided that the transition is 
sharp enough. The existence of such states in the regions 
of space that are in the midst of a change of topological 
charges has been extensively studied and form part of the 
physics of many systems. 

At the heart of the matter lies the fact that, under the 
effect of the transverse electric field, the pseudo-spinor 
eigenfunction warps the Brillouin zone in a non-trivial 
way enclosing it an integer number of times. This num- 
ber is a topological charge (the so-called Chern number) 
that cannot be untied in a continous fashion. Corre- 
spondingly changes in this topological charge are bound 
to be discontinuous and zero energy states are located 
nearby such discontinuities. The topological protection of 
the states together with their different effects from afore- 



mentioned effect has been suggested as a plausible valley- 
filter device. 

In the present work we include in the behavior of the 
chiral bands the effects of trigonal warping. We show 
that such distortion of the electronic bands in bilayer 
graphene dramatically modifies the behavior of the one- 
dimensional confined modes due to an inhomogeneous 
bias that changes sign across a channel. We have found 
up to four zero modes per valley per spin depending on 
the angle and the width of the domain wall. The effect 
is enhanced in wide domain walls however in such a case 
more branches can be found inside the gap. These other 
branches are localized within the domain wall and are of 
the same nature that the ones obtained in modelling a 
contact of two materials with inverted bands^.They are 
present regardless of the trigonal warping. 

A low energy reduced (dimensionless) Hamiltonian for 
the stacked AB bilayer, including the interaction respon- 
sible for the trigonal warping and an external field can 
be written as^ii^: 
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where Eq = (73/70)^71 ~ 4 meV, the units of k are fco = 
27371 /("\/3a7o), where we have used 70 = 3.16 eV, 71 = 
0.39 eV, and 73 0.315 eV, and a = 0.246 nm is the 
lattice constant. Additionally we employed the artificial 
parameter s to include continuously the effect of trigonal 
warping from s = 0, where it is strictly absent, to s = 1 
where it acquires its full value. The symbol k± stands for 
^kj. ± iky, where ^ = 1 for T-l around K' and ^ = — 1 for 
K. We can also write the Hamiltonian in the following 
manner ready to assess topological invariantsii 



H = £0 S{kx,ky) ■ a, 



where g{kx,ky) 
tion go(fc:E, ky) -- 



I {kx 7 ky 



y,V) and g^(fc^,n,y 



(2) 

In this equa- 
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{—kx,ky,0) where a wee the Pauli matrices. g,„ is the 
part responsible for the trigonal warping. 
The topological invariant N3 which describes K points 
is the winding number of the mapping of the sphere (T2 
around the K points to the 2-sphere of the unit vector 
g = g/|g| isii: 
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The topological charge can be found with no difficul- 
ties from the equations given above. It turns out that 
= £, sign(l/), independently of s. Despite this invari- 
ance the effects of trigonal warping are far from trivial. 
The trigonal warping splits the Fermi point, K into four 
pockets: three Dirac points away from the K point with 
charge -|-^ and the K point with charge — 
A reversal of the bias voltage over a region of space is 
therefore associated with a change of {SN^l = 2 for each 
valley and spin. This gives the minimum number of zero- 
energy left movers minus right moversii, there should be 
then at least two branches per valley per spin that cross 
the Fermi-level. 




FIG. 1. A transverse bias voltage with a domain wall that 
makes an angle i9 with the V — K axis. Propagative states 
are localized around the domain wall. In the figure the cir- 
cular spots represent carbon atoms arranged in a lattice with 
trigonal symmetry as in BLG. The hexagon represents the 
first Brillouin zone of BLG with the F point centered. The 
domain wall, between transverse potential V and —V has its 
geometry referred to the F — K axis, additionally, its width is 
denoted as £. 

We next proceed to model the electronic states con- 
fined by a bias-domain wall oriented along an axis h 
that makes an angle -d from the F — K axis. For that 
purpose we consider a potential of the form e()V{x) ~ 
EoVq tanh(a;j^/€) where x± is the coordinate perpendicu- 
lar to the h axis,^ is given in a (lattice constant) units, 
Vq is an adimensional measure of the strength of the po- 
tential and £0 is defined below Eq. ([T|). The domain 
wall therefore connects two regions of space with different 
topological characters. The change in the Chern number 



as one crosses the wall is 2. According to the index the- 
orem, this gives the number of zero-energy left movers 
minus right movers^i, there are at least two branches per 
valley per spin that crosses the Fermi-level. In order to 
study the detailed structure of those zero energy bands, 
we will in the following diagonalize the Hamiltonian. The 
diagonalizatiou of the Hamiltonian can be accomplished 
by noting that the momentum along the axis defined by 
the domain wall, h, is a good quantum number and states 
can be labeled by it, py. We then rotate the axis by an 
angle d such that the j/-axis lies along the domain wall. 
The Schrodinger equation can be transformed into: 



(4) 



where = i^(j)ya^ + l/2 se'^^^'^-) and 6 = 1/4 s^e'^^^'^- . 
The special case of s = (neglecting the effects of trigo- 
nal warping) can be solved analytically in terms of hyper- 
geometric functions)^ However, the inclusion of trigonal 
warping makes it more convenient to pursue a numerical 
approach. 

The resulting eigenproblem can be solved indepen- 
dently for each value of py. For a given value of py 
we discretize the axes across the wall and perform a 
numerical diagonalizatiou of the resulting Hamiltonian. 
The results are listed in Fig. © for several values of 
the angle -d {Va — 0.5). It is evident that despite the 
severe deformation suffered by the bands, there are still 
two bands that cross the Fermi level in each case. The 
velocity is heavily distorted and the trigonal distortion 
induces local energy extrema that are not present in the 
isotropic case. The location and depth of those valleys 
formed by trigonal distortion depend strongly on the 
angle that the wall makes with the lattice. 
We found that for an angle of 30° (90° for the other 
branch) the distortion is maximum, and is more evident 
from Fig. ^ that one of the confined bands crosses three 
times the zero energy. There is an angle (^ 23°(37°)) 
where two crossings also can be found making a total of 
three zero modes. The effect is enhanced by increasing 
the domain wall width. This interesting regularity can 
be found from the picture of fractionalization of Chern 
numbers. The integrand in eq. |3] (Berry Curvature) 
depends on the applied bias, the result is always the 
same but four well defined peaks can be found in the 
limit of small bias, one central dip with a contribution of 

— 1/2 and three others with a -1-1/2 value each, farther 
from the center and narrower the lower the bias. The 
total value of the integral is still +1 and on the other 
side of the wall, the values are opposite. 

The index theorem sets an inferior bound to the num- 
ber of zero modes in (iVg — iV|). If we sum the sign of the 
slopes at the crossings we obtain the number predicted 
by the theorem. Rotation of the domain wall breaks the 
symmetry between the two regions in such a way that 
the effective Chern numbers (A^3 ''') take different values 
on both sides of the wall. For instance (3 x (-1-1/2) — 1/2) 

- (3 X (—1/2) + 1/2) gives two, but the number three can 
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be obtained with (3 x (+1/2) + 1/2 - 1 /2) - (3 x (- 1 /2)) 
or shifting the ±1/2 to the other side, and the four zero 
modes can be formed with (3 x (+1/2) + 1/2) — (3 x 
(—1/2) — 1/2). Such reorganization of the topological 
charge on each side of the wall lead us to the proper 
number of zero modes. It is not possible to do a combi- 
nation to obtain more than four zero modes with those 
fractions. These half integer contribution has been ex- 
plored theoretically recently in several model a^'^i^^ . 
There are three contributions to the appearance of these 
zero modes, at low bias, important to split the topolog- 
ical charge, second an angle of the domain wall around 
30°(90°), breaking the symmetry, and a wide domain wall 
with a wide wave function and in turn a k vector more 
localized, the latter needed to take account of the slight 
splitting of the K point. 

Another important effect that can be appreciated is the 
lack of electron-hole symmetry in the system except in 
the case -d = 0° 60°.) This effect is expected since the 
Hamiltonian in Eq.Q fails to commute with the charge 
conjugation operator. 
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FIG. 2. One dimensional channels band structure for differ- 
ent angles of the domain wall. Trigonal distortion is seen to 
induce several effects on the band structure. All the plots are 
referred to the K' valley. ( Vb = 0.5 and 1 = 1). In (a) the 
bands for s = are included. Energies are given in eo units. 

The propagating modes need to lie in the channel 
defined by the domain wall, otherwise their energy will 
be gapped. In Fig. 2] it is clear that the singularities 
associated with the distortion of the confined channels 
lie within the gap and they are well localized inside 
the domain wall, singularities are different for different 
angles as expected from the band structure, we expect an 
increase in the conductivity with respect to the isotropic 
casein in the ballistic regime, due to the appearance of 
the extra zero mode. 

We have studied the distortion that trigonal warping 
exerts over the topologically confined modes in trans- 
verse bias-domain walls in BLG. We have shown that. 




FIG. 3. Band structures for different widths of the domain 
wall and for an angle of 30°. The four zero modes are evident 
and the other branches inside the gap in (b). Energies and £ 
are given in eo and a units respectively 
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FIG. 4. Local density of states, A{E,x), depicted as with a 
color map (arbitrary units) as function of e (vertical axis in eo 
units) and x (horizontal axis in a units)(Vo = 0.5 and I = 1). 
In (a) = 0° and (b) i9 = 30° . It is evident from these plots 
that the propagating states that cross the Fermi energy are 
located within the domain wall. 

despite their topological robustness, the band structure 
of those modes are heavily altered by the anisotropy in- 
duced by trigonal warping; depending on the angle and 
on the width of the domain wall up to four zero modes can 
be found, they can be explained in terms of a fractional- 
ization and reorganization of the topological charge. 

Acknowledgments ESM thanks L. Brey and L. Chico 
for helpful discussion. ASN work was partially funded 
by Proyecto Fondecyt 11070271 Proyecto Basal FB0807- 
CEDENNA and NUCLEO MILENIO N°P10-061-F. 
ESM acknowledges UTFSM for the internal grant PIIC- 
DGIP and Proyecto Basal FB0807-CEDENNA. PV ac- 
knowledges support from FONDECYT grant 1100508 
and CEDENNA, Chile. 

^Di Xiao, Wang Yao, and Qian Niu. Physical Review Letters, 

99(23):l-4, December 2007. 
^Henning Schomerus. Physical Review B, 82(16):l-6, October 

2010. 

■^Edward McCann and Vladimir I. Fal'ko. Phys. Rev. Lett., 

96(8):086805, Mar 2006. 
*Ivar Martin, Ya. M. Blanter, and A. F. Morpurgo. Phys. Rev. 

Lett., 100(3):036804, Jan 2008. 
^L. J. P. Xavier, J. M. Pereira, Jr., A. Chaves, G. A. Farias, and 

F. M. Peeters. Appl. Phys. Lett., 96(21):212108, 2010. 
6 J. C. Martinez, M. B. A. Jalil, and S. G. Tan. Appl. Phys. Lett., 

95(21):213106, 2009. 



4 



''Matthew Killi, Tzu-Chieh Wei, Ian Affleck, and Arun 
Paramekanti. Phys. Rev. Lett., 104(21):216406, May 2010. 

*J. R. Williams, L. DiCarlo, and C. M. Marcus. Science, 
317(5838) :638-641, 2007. 

9B. a. Volkov and O. A. Pankratov. JETP Lett, 42:178-181, 1985. 
lOMikito Koshino and Tsuneya Ando. Phys. Rev. B, 73{24):245403, 
Jun 2006. 

^^G. E. Volovik. The universe in a helium droplet. Oxford Univer- 
sity Press, Oxford, U.K. 2003. 



J. L. Maries, F. Guinea, and Maria A. H. Vozmediano. Phys. 
Rev. B, 75(15):155424, Apr 2007. 

Haruki Watanabe, Yasuhiro Hatsugai, and Hideo Aoki. Phys. 
Rev. B, 82(24):241403, Dec 2010. 

Takao Morinari. Phys. Rev. B, 72(10): 104502, Sep 2005. 
Michael Barbier, Panagiotis Vasilopoulos, and Frangois M. 
Peeters. Phil. Trans. R. Soc. A: Mathematical, Physical and 
Engineering Sciences, 368(1932):5499-5524, 2010. 



